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Abstract. We present a calculation of exciton states in semiconductor coupled 
quantum wells (CQWs) in the presence of electric and magnetic fields applied 
perpendicular to the QW plane. The exciton Schrodinger equation is solved in real 
space in three dimensions to obtain the Landau levels of both direct and indirect 
excitons. Calculation of the exciton energy levels and oscillator strengths enables 
mapping of the electric and magnetic field dependence of the exciton absorption 
spectrum. For the ground state of the system, we evaluate the Bohr radius, optical 
lifetime, binding energy and dipole moment. The exciton mass renormalization due 
to the magnetic field is calculated using a perturbative approach. We predict a non- 
monotonous dependence of the exciton ground state effective mass on magnetic field. 
Such a trend is explained in a classical picture, in terms of the ground state tending 
from an indirect to a direct exciton with increasing magnetic field. 
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1. Introduction 

Spatially-indirect excitons in coupled quantum wells (CQWs) present a model system 
for the study of a statistically degenerate Bose gas in solid-state materials. This is due to 
the long optical lifetime of indirect excitons in comparison to their thermalization time 
which permits the creation of a cold and dense exciton gas PEIE]. The exciton lifetime 
depends on the overlap of electron and hole wave functions. For indirect excitons, this 
can be varied over several orders of magnitude HE] by changing the electric held applied 
perpendicular to the growth direction. 

Numerous studies on indirect excitons have focused on Bose-Einstein condensation 
P d El E], spatial pattern formations uni EH Ea and spin transport phenomena 
[13]. Transport of indirect excitons has been studied extensively. Their dipolar nature 
provides a probe of the density via a blue shift in the emission [m ESI IB]- 
facilitates transport of excitons over tens of micrometers ini as indirect excitons can 
screen the QW disorder potential due to interface roughness and defects [13] . Moreover, 
the dipole-orientation of the excitons provides the possibility to control transport via 
patterned electrodes positioned adjacent to the CQWs OEniin]. Exciton transport 
has been studied in electrostatic traps [22l EH EH ESI EEji linear potential gradients 
[ini El] and stationary [27] and moving [28] lattices. The high degree of control led 
to demonstrations of devices such as excitonic optical transistors m- The properties 
of indirect excitons are also important for understanding a new breed of quasi-particles 
known as dipolaritons [30] • Such states are realized by embedding CQWs inside a planar 
Bragg-mirror microcavity [31] • At the resonant tunneling condition, the dipolariton 
is a superposition of a cavity photon, a spatially-direct and an indirect exciton and 
thus acquires a static dipole moment [321 EH- Dipolaritons have been suggested for 
observation of super radiant THz emission [33] and continuous THz lasing [35] . 

The application of magnetic helds to the CQW structure provides an extra degree 
of control. In exciton systems, magnetic helds are used for trapping [331 EH l38] - 
control of superhuidity [33 HD], tuning of exciton resonances [H] and studying spin- 
splitting phenomena [321 EH- The focus of this paper is the exciton effective mass 
enhancement due to a magnetic held [33]- From a qualitative view point, this ehect can 
be understood by considering a simple classical analogy. The Lorentz force experienced 
by two oppositely charged particles that travel perpendicular to a uniform magnetic 
held acts to increase their relative separation. This comes at an energy cost as it is 
opposed by the Coulomb attraction between the charges. The energy cost manifests 
itself as a resistance to change in momentum and thus an ehective mass increase, 
dependent on the applied magnetic held strength. We consider the case where both 
electric and magnetic helds are normal to the QW plane. By hnding a numerically 
exact solution of the exciton Schrodinger equation, we perform a precise calculation 
of the mass enhancement. This description is an important ingredient for simulating 
exciton transport kinetics in a magnetic held where the exciton ehective mass is a critical 

parameter [21 EH EH EH EH E5]. 
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Several works (Ha HTj for example) used variational methods to study indirect 
excitons in magnetic fields. However, these all lacked a calculation of the effective mass 
enhancement. A theoretical study of excitons in CQWs with hnite barrier and well 
width was done in [38], using two-dimensional (2D) free electron-hole pair states in 
magnetic held as a basis for expansion of the exciton wave function. This approach is 
only well suited for treating high magnetic helds. Also, a calculation of the exciton mass 
renormalization was absent. The mass renormalization in the 2D limit was calculated 
for excitons in single QWs Ha and indirect excitons in double QWs [50]. In those works, 
the electron and hole wave functions in the QW growth direction were approximated by 
delta functions located at the QW positions. This corresponds to the limit of inhnitely 
deep QWs of zero width and is applicable only for the analysis of either purely direct or 
purely indirect excitons in the case of high electric held and strong QW conhnement. It 
cannot describe in any detail the smooth crossover between direct and indirect exciton 
states at intermediate electric helds. 

The recently developed rigorous multi-sub-level approach [5] (MSLA) was 
previously used to calculate the CQW exciton states in perpendicular electric helds. The 
outcomes, which included a calculation of the dependence of the absorption spectrum 
and optical lifetime on electric held, were in good agreement with available experimental 
data. Here, we use the MSLA to solve the exciton Schrodinger equation in CQWs with 
both electric and magnetic helds applied perpendicular to the QW plane. This method 
solves the exciton in-plane problem by a numerical discretization in real space and is 
thus equally well suited for both low and high magnetic helds. We have rehned the 
MSLA to a highly accurate scheme by using Numerov’s algorithm [51], a fourth-order 
linear multistep method, for obtaining the energies and wave functions of the exciton 
states. It is found that the direct-to-indirect exciton crossover with varying electric 
and magnetic helds is an important ingredient in describing the magnetic held induced 
exciton ehective mass renormalization, which was not captured in earlier works. We note 
that the MSLA can also be used to solve the exciton Schrodinger equation coupled to 
Maxwell’s equations in order to model microcavity polaritons [33] • However, modelling 
polaritons in magnetic helds is beyond the scope of this paper and forms the subject of 
future study. 

In Section we describe the Schrodinger equation for a CQW exciton in the 
presence of electric and magnetic helds pointing in the growth direction and provide 
details of the MSLA. The perturbative approach used to obtain the electric and magnetic 
held dependence of the exciton ehective mass is outlined in Section In Section 
we present calculations of exciton states and their associated properties for 8-4-8 nm 
GaAs/AlGaAs CQWs. Comparisons with 2D methods are made in Section and we 
show that our method converges to these in the limit of zero QW thickness. Conclusions 
are made in Section 1 
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To calculate the bound exciton states in a semiconductor nanostructure with external 
bias and magnetic field orientated normal to the QW plane, we consider the following 
Hamiltonian: 

H{re, Th) = He{re) + Hh{rh) + Vc(|re - r/i|) + Eg, (1) 

with 

He,h{r) = Pe,h{r)^mfl{z)pe,h{r) + Ue,h{z), ( 2 ) 

Pe,ft(r) =-Vr ±-A(r), (3) 

t c 

Vc{r) = . (4) 

er 

Here, are the electron and hole coordinates. For a multi-layered heterostructure, 
we adopt a cylindrical coordinate system (p, z) with the z-axis along the QW growth 
direction. The electron and hole effective mass tensors, fh(,,h{z), are composed of in¬ 
plane and perpendicular components, and mfj^{z), respectively. As in [5], we 

neglect here any z-dependence in the in-plane electron and hole effective masses, which 
is justified by relatively low mass contrast in the heterostructures treated here and 
a minor contribution to the exciton problem of the electron and hole wave functions 
outside the well regions. me,h{z) are layer-dependent step functions, as are the QW 
confinement potentials {z). The confinement potential and the external bias are 
included in the potentials Ue,h{z) = {z) ± eFz where F is the electric field. For 

the band gap Eg and in-plane masses values for the QW layers are used, e is 

the average permittivity of the structure. For a perpendicular magnetic field (along the 
growth direction), one takes the vector potential, using the symmetric gauge, in the 
form A(r) = x r. 

It can be shown |S21 153] that the solution to the equation = ET has the 
following variable-separable form: 

^(re,rfe) = Xp(p,R)e*™'^(p(p,Ze,^ft) (5) 


with the center of mass motion of an exciton carrying 2D momentum P described solely 
by 

R 

h 


Xp(p, R) = exp <j i P -F -A(p) 


( 6 ) 


where R = p^me/M^ + pjpn'l/Mx and p = p^, — Pi^ = (p, 0) are the 2D in-plane center of 
mass and relative coordinates, respectively, Mx = me + is the in-plane exciton mass, 
and m is the exciton magnetic quantum number. Such a substitution enables removal 
of the center of mass coordinate from the problem. The exciton relative motion is then 
described by the Hamiltonian Hf derived in Appendix A For zero in-plane momentum. 
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Figure 1. Ground (el and hi) and first excited (e2 and h2) electron and hole states 
for (a) F = 3kV/cmand (b) F = 24kV/cm in 8-4-8-nm GaAs/Alo.asGao.erAs GQWs. 


the Hamiltonian ° takes the form 

z,,p) = + Hie,) + K(p) + n,(p) 

(\/p^ + 

where the Hamiltonians of the electron and hole perpendicular motion are 

d 1 d 


H. 


e,h 


.z = 


+ Ue^h{z). 


(7) 


( 8 ) 


2 dzmff^{z) dz 

These describe the Coulomb-uncorrelated single particle electron and hole states in the 
absence of a magnetic held, with the wave functions where the index q labels the 

electron (hole) states quantized in the CQW heterostructure potentials. The Schrodinger 
equations are solved numerically using the shooting method with 

Numerov’s algorithm. Figure shows the ground and hrst excited electron and hole 
states for hnite electric helds in a GaAs/AlGaAs GQW. 

The kinetic term in ([^, K{p), is given by 


r 


m 


(9) 


1_S_ 

p dp 

where p is the in-plane reduced exciton mass, 1/p = 1/ml + The potential due 

to the magnetic held is 

ehmB e^B'^p^ 


Vb(p) = 


2xc 


8pc'^ 


( 10 ) 


where 1/x = 1/mI — l/mf is the inverse mass diherence of the electron and hole, called 
the exciton magnetic dipole mass 
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The Schrodinger equation H^\k,m) = Ek,m\k.,m), in which \k,m) is an exciton 
state in a CQW in hnite electric and magnetic helds, is solved using the MSLA developed 
in |5]. The principle of the method is that we expand the exciton wave function 
gD(p,Ze,Zh) in ([^ into the set of Coulomb uncorrelated electron-hole pair states leading 
to 

|A^,m) = e‘”* 5^ (11) 

n=l 


where <hn(^e, Zh) = ipp^{ze)ip^^{zh)■ The mapping n i—)■ (p„, leads to Nf>Nh pair states 
formed from iVg electron and Nh hole states. We then calculate the radial components 
of the wave function using a matrix generalization of the shooting method 

with Numerov’s method incorporated in the hnite difference scheme. This solved the 
problem of the electron and hole in-plane and perpendicular motion mixed by the 
Coulomb interaction, leading to the exciton quantization labeled here by the index 
k. The numerical solution is generated on a logarithmic grid in the range p G [po,R]- 
A sufficiently small inner limit po is chosen so that the wave functions and extracted 
data have converged with respect to decreasing po. The outer limit R of the wave 
functions is sufficiently large that all radial components asymptotically approach zero. 
The existence of such a limit is guaranteed for non-zero B as the potential Vb grows 
with p^. This is in contrast to the case of i? = 0 where truncation of the solution range 
leads to a discretization of the unbound e-h continuum states. For the range of electric 
helds considered {F < 25kV/cm), it was previously found for B = 0 that high accuracy 
is achieved using two electron and hole states (A^"e = = 2) [5]. 

The state with zero angular momentum has radiative linewidth T^^ = 

{7ie‘^hf^’‘^)/{y/emoc). The oscillator strength, is calculated from the overlap of 
the electron and hole wave functions. 


/“> = 


2moEkp\dc 

h? 






dz 


n=l 


( 12 ) 


where dev is the dipole matrix element between the conduction and valence bands. 


3. Exciton mass renormalization in a magnetic field 


For an exciton with a hnite center of mass momentum P 7 ^ 0, the Hamiltonian Q is 
modihed to 


Hf = m + Fp 


where the perturbation Vp is given by (see Appendix A) 

p2 


Fp = 


2 e 


2iVA AAc 


P ■ A(p), 


(13) 


(14) 


and P = |P|. To calculate the exciton mass renormalization in a magnetic held, we 
treat P as a small parameter and use perturbation theory up to 2nd order. The 1st 
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term of (14) contributes to the exciton energy in 1st order only leading to the bare mass. 
The 2nd term, vanishing in 1st order, gives rise to a quadratic in P correction in 2nd 
order and is thus responsible for the mass renormalization. Neglecting non-parabolicity 
of the exciton band, which would be accounted for by treating the 2nd term in higher 
perturbation orders, we find the correction to the exciton energy proportional to and 
renormalized effective mass of exciton state \k,m) in magnetic field: 




k,m 


1 


2e 


M^cP 


E 

s = ±l 


|(fc, m |P-A(p)| j, m + s)|' 

^k,m 


(15) 


Here the index j counts over the eigenstates with the angular momenta mil which 
only contribute in 2nd order and which are calculated at given values of electric and 
magnetic held. Details of the calculation of the matrix elements in (15) are given in 
Appendix B 


The advantage of our approach compared to some previous calculations of the 
exciton mass renormalization is that the perturbation is in P only, allowing the H-held 
to be arbitrarily large. In contrast, the approach developed in |1H] used the magnetic 
held as a small parameter of the perturbation theory. Consequently, the applicability 
of the latter is restricted to low magnetic held (up to ~ 2 T in the structure considered 
in this paper, see also the comparison in Section]^. Another signihcant beneht of our 
approach is that the full 3D solution of the exciton Schrodinger equation ([^ describes 
the inter-well coupling that is neglected by any 2D approximation [4U1130] . 


4. Application to GaAs/AlGaAs CQWs 

All results in this section refer to the structures studied in [D El n [HI EQ EH EB] 
which consist of symmetric GaAs/Alo.ssGao.eTAs QWs with barrier and well thickness 
of 4nm and 8nm, respectively. A complete list of parameters used in the calculations 
are given in table [T} Although we restrict the present analysis to two QWs, the method 
is applicable to any number of QWs since one can always obtain, either analytically or 
numerically, eigenstates of the single particle electron and hole Hamiltonians in ([^. 

Figure a) shows the optical transition energy for different exciton states as a 
function of electric field at i? = 10 T. For each state, the circle area is proportional to 
the oscillator strength. One can identify four different exciton states with each state 
being duplicated in each of the first three Landau levels. Firstly, we see a pair of bright 
states at E^ — Eg ^ 51meV whose energy is almost independent of electric field. These 
are the first Landau level (Li) direct exciton states occupying adjacent QWs. The 
direct exciton states are close in energy and appear almost as a single state in figure 
l^a). Secondly, we see a pair of darker states, one whose energy decreases and one 
whose energy increases with electric field from E^ — Eg 60meV at F = 0. These 
are the Li indirect exciton states. Their transition energy dependence on electric field 
is almost linear with the gradient corresponding to the nominal distance between the 
QWs. The state with energy that decreases linearly with F is composed predominantly 
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Table 1. Parameters of the model 


£ 

Relative permittivity 

12.5 

E, 

GaAs band gap 

1.519eV 

mf{z) 

Electron mass in QW 

0.0665 mo 


Electron mass in barrier 

0.0941 mo 

mf{z) 

Hole mass in QW 

0.34 mo 


Hole mass in barrier 

0.48 mo 

M, 

In-plane exciton mass 

0.22 mo 

p 

In-plane reduced exciton mass 

0.042 mo 

K 

Exciton magnetic dipole mass 

0.15 mo 

^cv 

Dipole matrix element 

0.6 nm 

Po 

Inner wave function boundary 

0.025 nm 

R 

Outer wave function boundary 

500 nm 


Number of in-plane grid points 

300 


z-gcid spacing 

0.1 nm 


of the ground electron and hole states elhl shown in hgure The other is composed 
predominantly of the excited electron and hole states e2h2. The corresponding L 2 and 
L 3 replicas of these four states appear at higher energies. The spacing between the 
Landau levels for each type of exciton is approximately equal. 

In the previous work using the MSLA |5], which did not include magnetic held, 
it was necessary to calculate a large number of unbound e-h pair states in order to 
determine the exciton spectrum. These states were discretized due to a hnite domain 
of the solution which imposed an inhnite potential barrier at p = R. It was found that 
with increasing R, the density of such discretized continuum states increased while their 
oscillator strengths decreased. The absorption spectrum converged for sufficiently large 
R. In the present case, however, the magnetic held removes any need to consider the 
continuum states. The potential due to the magnetic held, (10), contains a term in 
which naturally closes the domain of the solutions. Therefore, a hnite value of R can 
always be found such that the wave functions are sufficiently small there. This point 
explains the apparent cleanliness of hgure [^a) which does not contain any unbound e-h 
pair states. 

Figure [^b) shows a magnihcation in the region of anti-crossing between the Li 
direct and indirect exciton states, labeled in hgure [^a). We identify the electric held at 
the anti-crossing Fa-c as the value of F where the minimum separation occurs between 
the upper and lower black curves. This is shown in hgure [^c) as a function of B for the 
hrst three Landau levels. The magnetic held provides a means to tune Fa-c over several 
kV/cm. The corresponding energy splitting between the upper and lower black curves 
in hgure [^b) is about 2 meV and varies by less than 10% in the investigated magnetic 
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Figure 2. (a) Electric field dependence of the optical transition energy for different 
exciton states. The circle area is proportional to the exciton oscillator strength, (b) 
Magnified image [indicated in (a) by the boxed area] of the anti-crossing between the 
direct and indirect exciton states in the first Landau level, Li. Black curves show 
the exact transition energies. The red arrow shows the minimum energy separation 
between the lower and upper curves, (c) The electric field, ik-c, at which the minimum 
separation occurs between the upper and lower curves in (b) as a function of B. 
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Figure 3. Magnetic field dependence of the exciton absorption spectrum at (a) 
QkV/cm and (b) 24kV/cm. The spectra are convoluted with a Gaussian of width 
ImeV. Both spectra are normalized to 1 at the maximum. 


field range. 

Using the Lorentzian model of absorbing oscillators, the absorption spectrnm is 
calculated from the transition energy and linewidth of each exciton state. The magnetic 
field dependence of the absorption spectrum is plotted in hgure For clarity, the 
spectrum has been broadened by convolution with a Gaussian of full width at half 
maximum of 1 meV. Most clearly visible is the Landau fan of the direct exciton states. 
These bright spectral lines are due to the two direct exciton states overlapped in energy. 
Fainter fans are seen for the indirect exciton states - one starting from about 20meV 
and another from about 80meV in hgure [^b). The lower (upper) fan corresponds to 
indirect excitons composed mainly of an electron and a hole in the ground (excited) 
state, elhl (e2h2). 

A number of properties of each eigenstate can be calculated from the exciton wave 
functions (pTj). In hgure we focus on the characteristics of the ground state of the 
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Figure 4. Electric and magnetic field dependence of (a) the exciton ground state 
in-plane Bohr radius, (b) dipole length, (c) radiative lifetime, and (d) binding energy. 


system. At zero magnetic field, the results presented here agree |S1] with those of 
[5]. Figure]^ a) shows the electric and magnetic field dependence of the in-plane Bohr 
radius, = \/ (p^). Increasing the electric field increases due to a reduction in the e- 
h Coulomb interaction which is caused by their increased separation in the ^ direction. 
Increasing the magnetic held has the reverse effect and shrinks the wave function in 
the QW plane. This shrinkage enhances the e-h interaction which, in turn, leads to 
a more likely localization of the electron within the same QW as the hole and thus 
the ground state becomes a direct exciton. This is seen in hgure |^b) where a sharp 
reduction in the dipole moment \{ze — Zh)\ takes place for increasing B. The value of 
B at which this transition occurs increases with F as a stronger conhnement of the in¬ 
plane wave function is needed to sufficiently enhance the e-h interaction and to induce 
the indirect-to-direct transition. In the indirect exciton regime at high electric field 
(F = 24kV/cm), the dipole moment corresponds to approximately the center-to-center 
distance of the QWs which is 12 nm. 

From the radiative linewidth, the radiative lifetime of the ground state exciton, 
tr = V(2r^^), is found. Figure l^c) shows tr as a function of electric and magnetic 
fields. Tr is inversely proportional to the overlap of the electron and hole wave functions 
(see hgure [^. Therefore, we hnd a decrease in tr with decreasing F or increasing 
B, consistent with the shrinking of the exciton radius shown in hgures |^a) and|^b). 
We note that in experiments [SS], contrary to our results, the radiative lifetime of 
QW excitons is found to increase with increasing magnetic held. This, however, can be 
explained by considering the thermal distribution of excitons. Only exciton states inside 
the light cone, i.e. those with in-plane momentum P < P^, can decay optically. The 
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Figure 5. Electric and magnetic field dependence of the exciton ground state effective 
mass calculated using the full MSLA. Inset: The same calculated within the classical 
model described in text, using the data in figures |^a) and|^b). 


momentum Pj corresponds to the intersection of the exciton and photon dispersions. 
The increase in the exciton effective mass (see hgure causes a reduction in which 
reduces the fraction of excitons coupled to light. This increases the effective lifetime of 
the exciton gas—the quantity which is measured in experiments via the decay rate of 
the photoluminescence intensity mi- 

Figure ^d) shows the binding energy Eb = E^'’ — Eq^ + heB/ifiiic) where Eq^ 
is the calculated exciton ground state energy in magnetic held, = E^P E'f is the 
energy of the unbound e-h pair state in zero magnetic held, and the last term takes into 
account the Landau quantization energy of the free pair in magnetic held. The decrease 
in the size of the exciton, caused by either increasing B or decreasing F, leads to an 
enhancement of the e-h interaction [56]. This causes a tighter binding of the exciton 
and, therefore, an increase in Eb- 

The ehective exciton mass, is calculated using second order perturbation 

theory (see Section]^. Figure]^ shows the electric and magnetic held dependence of the 
ehective mass, Mq g, for the exciton ground state. We identify two limiting cases: Firstly, 
for small electric helds {F < 3 kV/cm) where the exciton ground state is predominantly 
direct, there is a monotonic increase of the ehective mass with increasing magnetic 
held. The rate of this increase is weakly dependent on the electric held. Secondly, for 
electric helds greater than 3 kV/cm, there is an initial increase in the ehective mass with 
increasing magnetic held, followed by an eventual decrease. At higher magnetic helds, 
the ehective mass asymptotically approaches the direct exciton ehective mass. 

To interpret the data presented in hgure we use the classical analogy introduced 
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in Section [T] The exciton is approximated as two oppositely charged masses at an 
equilibrium separation = \/r‘^ + {ze — ZhY with center of mass motion perpendicular 
to the magnetic held. The Lorentz force Fl = eBv/c acts on each particle to separate 
the charges but is balanced by the Coulomb restoring force Fc that is approximately 
linear for small changes in separation (corresponding to small velocity and small Lorentz 
force): Fl = Fc ~ knAp. In this picture, one derives from the energy dependence 
on the exciton velocity v, M^v‘^/2 + k^{Apfi/2, an estimate of the effective mass, 
= Mx + e‘^B‘^/{c^kji). The ‘spring constant’ k^ is the coefficient of the restoring 
force. It can be approximated via linearization of the slope in the Coulomb potential at 
p = tb which gives 


kn — — 


f)2 _ ■ 


(16) 


p=rB 


The inset in hgure shows this classical calculation using the parameters tb and 
{ze — Zh^ characterizing the exciton taken from the full calculation, see hgures |^a) 
and|^b). The curves for F = 3kV/cm and 6kV fully describe the main trends seen in 
the full calculation and even demonstrate to some extent a quantitative agreement. For 
larger electric and magnetic helds this classical model fails because the spring constant 
kn becomes too small and the linearization of the restoring force used in the model no 
longer works. 

The classical model allows us to understand the exciton mass dependence in hgure 
Indeed, from the calculation presented in hgure we hnd a greater k^ for direct 
excitons than indirect excitons. Qualitatively, this is because the indirect exciton has 
a weaker binding due to the increased spatial separation of the electron and hole. The 
diherence in leads to a greater rate of ehective mass increase with magnetic held for 
indirect excitons than for direct excitons. This explains the ehective mass dependence 
in the low and high electric held limits in hgure At intermediate helds (6 —12kV/cm), 
the initial increase and then decrease of the ehective mass is a consequence of the ground 
state tending from an indirect to a direct exciton. This happens due to a shrinkage of the 
exciton wave function caused by the magnetic held. In turn, the e-h Coulomb attraction 
is enhanced, making the direct state lower in energy than the indirect state. 


5. Comparison with the 2D limit 

In order to verify our model, we make comparisons with previous results in which the 
limiting case of 2D excitons were treated. To make this comparison, we include just 
one electron and one hole state [N^ = Nh = 1 in ( [IT| )] taking the wave functions 
in the form \'ip^{z)\‘^ = 6{z) and = S{z — dfi, where 6{z) is the Dirac delta 

function and dz is the nominal distance between the QWs. Technically, this is done 
using normalized Gaussian functions for 'if>‘^'^{ze,h), for the electron and hole centered on 
0 and dz, respectively, and taking the limit of their widths to zero. 

First, we consider the mass renormalization of 2D direct excitons in single QWs 
{dz = 0). This was calculated by Arseev and Dzyubenko [IH]. The analytical expression 
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for the exciton mass, which was derived by treating the magnetic held as a perturbation, 
was given as |57| 


M. 


2D 


1 


42/i 

162M„ 


(17) 


This is relevant in the low magnetic held limit where the magnetic length is much greater 
than the exciton Bohr radius, Ib = \/hc/eB S> Ox = efiffue^. For the considered 
system, this corresponds to i? -C 2.6 T. Figure a) shows the inverse of the ehective 
mass renormalization as a function of calculated numerically via (15) and using the 
analytic expression (0. There is a strong agreement between the diherent approaches 
for small B. For increasing B above IT, the results depart as expected, since the 
magnetic length approaches the Bohr radius. 

We also compare our results with the work of Lozovik [50] which treated the case of 
2D indirect excitons in CQWs. In that approach, the exciton Schrodinger equation was 
solved with a hnite center of mass momentum that enabled calculation of the dispersion 
relation. The effective mass was then extracted from a parabolic £t to the bottom of the 
dispersion band. The mass calculated by this method is shown in figure [^b). To make 
a suitable comparison, we calculate the 2D effective mass using the same parameters 
as in jSn] and [Hj: d^ = 11.5nm, /r = 0.049mo, x = 0.11 mo, and e = 12.1. The full 
calculation using the same parameters with F = 24kV/cm is shown for reference. 

In hgure [^b), we also show the effective mass dependence on magnetic held as 
was determined experimentally [201 SI] using an in-plane component in the magnetic 
held to map the indirect exciton dispersion by shifting it with respect to the photon 
cone. We note that the quantitative agreement between theory and experiment could be 
improved by further adjustment of the parameters /i and x. The calculated mass is quite 
sensitive to these parameters which are not known to great precision. Other possible 
sources of discrepancy could be the inhomogeneous permittivity, non-parabolicities of 
the bands and valence band mixing. One should also consider the limited accuracy 
of the effective mass measurement which entails htting a parabola to an almost hat 
dispersion. The dispersion curve is extracted from the energy of the exciton line in the 
emission spectrum. This can be inhuenced by the disorder intrinsic to the CQW sample 
and the density dependent blue shift due to the exciton-exciton interaction. On that 
basis, our data is in good agreement with the experimentally measured exciton ehective 
mass and the combination of the ehects that we neglect is rather minor. 


6. Conclusions 

We have studied the combined ehect of perpendicular electric and magnetic helds on 
exciton states in coupled quantum wells. This was done using the rigorous multi-sub- 
level approach (MSLA) to solve the full exciton Schrodinger equation in the ehective 
mass approximation. In this approach, the exciton wave function is expanded into 
Coulomb uncorrelated electron-hole pair states and found using a matrix generalization 
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Figure 6. (a) Present numerical calculation (solid red line) and analytic calculation 
of [in] using ( [TtI (dashed black line) of the inverse of the mass renormalization against 
for 2D excitons in a single QW. The inset shows the corresponding effective mass 
increase as a function of B. (b) Indirect exciton effective mass as a function of -B, as 
determined by the full calculation (green dashed line), using 2D electron and hole wave 
functions separated by dz = 11.5 nm (blue line), via the dispersion relation calculated 
in [SOI (black circles with a line) and experimentally jUj (red triangles with error bars 
and a line). The error bars show the errors extracted from the data presented in [44] . 


of the shooting method with high stability and accuracy provided by Numerov’s 
algorithm. We calculated the optical transition energy and linewidths for different 
exciton states and obtained the electric and magnetic held dependence of the exciton 
absorption spectrum. For the exciton ground state, we evaluated the in-plane Bohr 
radius, static dipole moment, radiative lifetime and binding energy. Furthermore, we 
have calculated the magnetic held-induced exciton mass renormalization in arbitrary 
electric and magnetic held, by treating the exciton momentum as a small perturbation. 
In particular, our calculations demonstrate a non-monotonous behavior of the exciton 
ehective mass in magnetic held which is understood in terms of the direct-to-indirect 
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exciton transition and also illnstrated in a simple classical model. Our results for the 
quantum well width tending to zero are in good agreement with different 2D approaches 
available in the literature. However, in contrast to the 2D approaches, the full MSLA 
captures the direct-to-indirect transitions of the exciton ground state for varying electric 
and magnetic helds which are necessary for a full and detailed description of the system. 
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Appendix A. Exciton Hamiltonian in a CQW in perpendicular electric and 
magnetic fields 


Using the factorizable form of the wave function we make a transformation of the 
full Hamiltonian ([^ 

(U, rh) = Ze, Zh)xp , (A.l) 

removing the in-plane center of mass coordinate from the problem. The exciton 
Hamiltonian then takes the form 


H^{p, z,, zi) = H^(z,) + Hi(zt) + WE(p) + E, 

+ Vc lx9‘ + W - Zhfi , 

where and Vq are given by ( 2 ) and 
Wl{p) = XF^HB{pe,Ph)XP 

with 

HBiPe.Ph) = ^{-Vl-2txA,-Ve + l^Al) 

2me 

+Aj (-Vj + 2»7A».Vi + . 

2ml 


(A. 2 ) 

(A.3) 


(A.4) 


For brevity of notation, we have introduced in (A.4) 7 = e/hc, Ae^ = A{p^ fi), and 
Ve 7 = where Vp is a 2D nabla-operator. For the same reason, we also drop 

below the index P in yp, x in M^, and || in Taking an arbitrary function f{p) we 

r 


obtain from (A.4): 


HBXf = 


2mf. 


xVf/-2VeX-Ve/-/Vfx 


2zxA,-ixVef + fVex) +x"Alxf 


2mh - 


-xVif - 2VhX ■ V,/ - fVix 
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+ 2i'yAh-{x^hf + f^hX) + 7 Axf 


Noting that Ve/ = V/ and Vh/ = —V/, (A.5) can be written as 


- 2V/ ^ 


m. 


ruh 


-/ (— + —) - 2i7— ■ (XV/ + /V.x) 

\ me mh J me 

+2i7— ■ (-\'V/ + /Vix) + 7 AA + 


rrih 


me 


mhj 


Finally, noting that 


V.7A(p) .R=tv.7B X 


VM' 


M 


Ph 


= -Ve,hB ■ PeX Ph = -Ah,e: 


obtain 




(A.5) 

(A.6) 

(A.7) 

(A.8) 


(A.9) 

(A.IO) 


and after some algebraic simplihcations using A(p) = Ag — A^, hnd the required term 


hFg of the exciton Hamiltonian (A.2): 


’h7p)--|v^-^A(p).VA^ 




(A.11) 


This term is the only one in (A.2) which carries in the dependence on the magnetic held 
B and the exciton in-plane momentum P. For the exciton wave function in the form of 
(|^, with the split-off angular motion, the exciton Hamiltonian is further reduced to ([^ 
for P = 0 and to (14) for P ^ 0. 


Appendix B. Perturbative calculation of the exciton mass enhancement due 
to a magnetic field 


We treat the exciton center of mass momentum P as a perturbation and hrst seek 
solutions in a hnite magnetic held with P = 0. Without loss of generality, we can 
take P along the y-axis, which gives P ■ A(p) = ^BPpcos{(f)) in polar coordinates, and 
evaluate the matrix elements. 


irff = -(k,m 

k,j P ' 


2 e 


APc 


P ■ A(p) 


m -|- s, j 


eB 

M^c 



r-2-K 


dzgdzh d(p p(ippcos(0)e“'^ 




n 


(B.l) 




















Exciton effective mass enhancement in electric and magnetic fields 


18 


We make use of the following identities: 


cos(0) dfj) = TT(5^2,1 

~ ^Pn,P„n 


' —OO 
POO 


'^q„ {^h) dZh ) 


dpn,p„'dqn,q„/ dn,n'- 


The matrix elements then simplify to 

Tm,s r TTcB 

-lu ^ = 0 q 2 1 - 




(B.2) 

(B.3) 

(B.4) 

(B.6) 

(B.6) 


The correction, to the effective mass of exciton state |fc, m) is then given by the 

2 nd perturbation order: 


2M, 


(B) 


E 


k,m 


(I, 


rn,-l\2 
k,j > 


E, 


+ 


(I, 


m,l\2 


k,m 


Ejm—l 


k,m 


E 


'j,ra-\rl 


(B.7) 


Finally, the renormalized effective mass of an exciton in a magnetic held is 

determined by 

■ ■ ■ 


Mt 




M, 


(B)' 


Figure |B^ shows some illustrative examples of the convergence of the effective mass 
with respect to jmax for A; = m = 0. Here, jmax is the number of states taken into account 


in the perturbation series, (B.7). It can be seen that the series converges rapidly for the 


range of electric and magnetic helds considered here, jmax = 15 was used in all data 

except hgure |Bl[ achieving a relative error in the effective mass of less than 10“'^. 
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